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Abstract 

This paper proposes to unify fading distributions by modeling the magnitude-squared of the in¬ 
stantaneous channel gain as an infinitely divisible random variable. A random variable is said to be 
infinitely divisible, if it can be written as a sum of n > 1 independent and identically distributed 
random variables, for each n. Infinitely divisible random variables have many interesting mathematical 
properties, which can be applied in the performance analysis of wireless systems. It is shown that the 
proposed unification subsumes several unifications of fading distributions previously proposed in the 
wireless communications literature. In fact, almost every distribution used to model multipath, shadowing 
and composite multipath/shadowing is shown to be included in the class of infinitely divisible random 
variables. 


1. Introduction 

Wireless eommunieations systems are primarily subjeet to two types of fading effeets. Small- 
seale multipath fading, corresponds to short-term variation of the signal, as it is the result 
of constructive and destructive combinations of random reflections and scattering of signal 
components. Models such as Rayleigh, Rician and Nakagami-m distributions have been proposed 
to capture the distribution of the envelope of the multipath fading random variable (RV). The 
second type of fading is called shadow fading, which captures the long-term variation of the 
mean signal level, and is typically modelled as a lognormal RV. In certain scenarios, a composite 
multipath/shadowing distribution is found to be appropriate, and is modeled as the product of a 
multipath fading RV and a shadowing RV. It is of interest to see if the typical distributions used 
for multipath, shadowing, and composite multipath/shadowing may be unified under a common 
class with desirable analytical properties. Through a unification of fading distributions, it may be 
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possible to obtain canonical expressions for the performance metrics of fading channels, thereby 
simplifying performance analysis. The unified model may also permit the comparison of two 
different fading distributions with respect to system performance metrics such as the average 
symbol error rate (SER), using stochastic orders, which are partial orders on RVs. Therefore, it 
is compelling to develop a unified study of fading distributions with these goals in mind. 

There have been multiple efforts toward unifying fading distributions in the past. Spherically 
invariant random processes have been proposed for the fading envelope as variance mixture of 
Rayleigh distributions |[I1 . An unification proposed in |l2l models the instantaneous channel power 
RV as the product of a gamma and a generalized-gamma RVs. While this model encompasses 
many known fading distributions, it does not include some distributions such as the Pareto 
distribution, which has been used in [l3]| to account for interference. Another unification in 
the literature is the generalized gamma distribution (also known as the Stacy distribution [|4|1. 
which is shown to correspond to the distribution of the envelope, when the received signal is 
composed of clusters of multipath components, and the receiver possesses power-law nonlinearity 
m. Reference [[6]| also proposes a fading model, however, it does not include some popularly 
used fading distributions such as the Rician distribution. For the first time in the literature, 
we propose to unify fading distributions by modeling their instantaneous channel power as an 
infinitely divisible (ID) RV. A RV is said to be infinitely divisible, if it can be written as a sum of 
n > 1 independent and identically distributed (i.i.d.) RVs, for each n. Intuitively, this means that 
the channel can be viewed as a (time/frequency/antenna) diversity combining system with n i.i.d. 
channels, for every n > 1. Infinitely divisible RVs have many interesting mathematical properties 
ifTl , which are relevant in the performance analysis of wireless systems. In this paper, almost 
every distribution used to model multipath, shadowing and composite multipath/shadowing will 
be shown to be ID RVs. Most of these distributions are members of a special subclass of ID RVs 
known as generalized gamma convolutions (GGCs). Consequently, special attention is devoted 
to this subclass of ID RVs. GGCs are defined as the limits of sums of Gamma RVs, potentially 
with different parameters, and have been thoroughly studied in IfSl- lfTTI . There are numerous 
applications of GGCs in financial economics and risk analysis ifTSl . [IT^ . GGCs possess 

remarkable closure properties under summations or products. These properties make GGCs an 
attractive model for multipath fading, shadowing or composite models that incorporate both. 
Despite possessing properties which are useful in fading system analysis, to the best of our 
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knowledge, no applications of GGCs has been found in the wireless communications literature. 

The class of GGCs includes a surprisingly large number of fading models such as the log¬ 
normal, Rayleigh, Nakagami-m, generalized gamma (Stacy), Pareto, beta, inverse Gaussian and 
the positive a-stable distributions. Moreover, composite multipath/shadowing models such as 
Rayleigh-lognormal, Nakagami-lognormal, and special cases of the spherically invariant random 
process, and the generalized model proposed in [[21 are also GGCs. From a wireless perspective, if 
a system has GGC distributed instantaneous channel power, then the fading channel is equivalent 
to a coherent linear combination of (possibly infinitely) many independent branches where 
the instantaneous fading power is gamma distributed. Unifying instantaneous channel power 
distributions as GGCs, it is possible to obtain novel generalized expressions for performance 
metrics of wireless systems, such as the ergodic capacity and the average SER. Furthermore, 
this unification reveals new stochastic ordering relations applicable to members of GGC. Lastly, 
the conditions under which systems such as diversity combining schemes have GGC end-to- 
end instantaneous channel power, are obtained. This facilitates one to obtain bounds on the 
performance of such systems. 

The rest of the paper is organized as follows. In Section HI] provides the required relevant 
mathematical background. The system model and the performance metrics under consideration 
are given in Section jllll A short exposition of infinite divisibility highlighting applications in 
wireless communications is presented in Section jlVl Section IIV-DI provides examples of fading 
distributions which are ID. Performance analysis of GGC fading channels is delegated to Section 
rvll Section jvnj describes the stochastic ordering of GGC RVs. In Section IVIIIl multi-antenna 
systems with maximum ratio combining (MRC), is considered. Finally, Monte-Carlo simulations 
are presented in Section IS 

A. Notation and Conventions 

The set of real numbers and positive integers are denoted by M and N respectively, while 
all other sets are denoted using script font. For a set S, the indicator function / [x G 5] = 1 if 
X G iS, and 0 otherwise. For a Radon measure /r(-) on M+, ^{u) is used to represent /i([0, m]). For 
a RV X, the cumulative distribution function (CDF) and the probability density function (PDF) 
are denoted by Fx (■) and /x (•) respectively. E [^((X)] is used to denote the expectation of the 
function g{-) over the PDF of X. The Laplace transform and Laplace exponent of a nonnegative 
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RV X are defined as (•) = E [exp(—sX)] and 7 x(-s) = — logE [exp(—sX)] respeetively, 
whieh are defined for s > 0. All logarithms are natural logarithms. We write fi{x) f 2 {x), as 
a; —)■ a to indieate Muix^a fiix)/ f 2 ix) = 1. 

II. Mathematical Preliminaries 

A. Bernstein Functions 

A funetion g : M’*' —)■ M’*' is said to be eompletely monotone (e.m.) if it has derivatives of all 
orders, whieh satisfy {—g{x)/(ix^ > 0, for k G MU{0}. Aeeording to Bernstein’s theorem, 
every e.m. funetion is characterized as a mixture of decaying exponential functions IflTl . Many 
performance metrics of communication systems such as SERs of M-QAM and M-PSK are e.m. 
functions of the average SNR ifTFI . 

If a function g : M+ —)■ is such that dg{x)/dx is e.m., then it is known as a Bernstein 

function [fTTIl . Every Bernstein function can be written in the form 

OO 

g{x) = a + bx + J ~ exp(—sx))r(ds) , (1) 

0 

for some a, 6 > 0, and a non-negative measure r on M+ which satisfies sr(ds) + r(ds) < 
OO. In this case, r is known as the Eevy measure of g. Bernstein functions are relevant to 
this work, because the Eaplace exponents of non-negative infinitely divisible RVs introduced in 
Section UYI are Bernstein functions |[l7l p. 37]. 

The subset of Bernstein functions, for which st(s) is e.m., constitutes the set of Thorin- 
Bemstein functions, TBX. Any function in this set admits the following representation, in 
addition to O: 

CXD 

g{x) = a + hx + j log(l -f sx)fi{ds) , (2) 

0 

for some non-negative a, b, and a non-negative measure on M+ which satisfies | log s|/i(ds)-|- 
sp(ds) < OO. In this case, /x is called the Thorin measure of g. In this work, we refer to 
/i(ds) as the mass of the Thorin measure, or the Thorin mass. It will be seen that the Thorin 
mass of the instantaneous channel power is the diversity order. We will henceforth refer to the 
Thorin mass as diversity order. 
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An example of a Thorin-Bernstein funetion relevant to wireless eommunieations is g{x) = 
log(l + a;). This is beeause, if pX is viewed as the instantaneous SNR of a ehannel, where p is 
the average SNR, then E [g{pX)] represents the ergodie eapaeity of the ehannel. 

B. Regular Variation 

Some eoneepts from regular variation theory will be used in Seetion IIV-B[ Intuitively, regular 
variation eaptures asymptotic polynomial-like behavior. 

Definition 1. A real valued function H{x) : M’*' M’*' is said to be regularly varying at oo (at t)) 

if Ywix^oo H(tx)/H(x) (Xvaix^a H(tx)/H(x)) exists, and is equal to f, with 0 < |r| < oo. It is 
said to be slowly varying at oo (or at D) if r = D. If H is differentiable with h(x) = dH(x)/dx, 
and h(x) = mx"^~^l(x), where I is slowly varying at 0, then H(x) ~ x"^l(x), as x ^ 0. 

C. Stochastic Ordering 

Stochastic orders are binary relations defined on probability distributions which capture intu¬ 
itive notions like being larger or being more variable. A systematic development of stochastic 
ordering theory can be found in [|T9ll . Some stochastic orders relevant to this paper are now 
presented. In what follows, we define stochastic orders on measures which might not have unit 
mass, and therefore are not necessarily probability measures, since this generality will be needed 
in Section IVIIi 

1) Convex Order: Let X and Y be two random variables such that 


E |9(A')] < E |g(F)] 


(3) 


for all convex functions g, provided the expectations exist. Then X is said to be smaller than 
Y in the convex order (denoted as X <cx Y). 

2) Laplace Transform Order: Let /i be a nonnegative measure. Its Laplace transform is defined 


as exp(—sM)/i(ds), for m > 0. This is used to define the LT order between two nonnegative 


measures. If pi and p 2 are two nonnegative measures, then pi is said to be dominated by p 2 in 
the Laplace transform sense (denoted by pi <Lt P 2 ) whenever: 



( 4 ) 
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This definition of the LT order is somewhat more general than definition found in [fT^ Chapter 
5], whieh applies to only probability measures. When /xi and are probability measures 
eorresponding to RVs X and Y respeetively, useful properties are obtained. E [^'(X)] > E [g{Y)] 
for all c.m. funetions g{-) is equivalent to X <Lt Y IfT^ pp. 96]. A similar result with a reversal 
in the inequality states that E [g{X)] < E [g{Y)] for all Bernstein funetions g{-) is equivalent to 

A <Lt 

3) Shannon Transform Order: In what follows, a generalized version of the ergodie eapaeity 
order [|^ applieable to order nonnegative measures is deseribed. 

Let /i be a non-negative measure on R+. Its Shannon transform is defined as Jg°°log(l + 
p/s)/r(ds). If /ii and fi2 are two non-negative measures for whieh the Shannon transforms exist 
Vp > 0 , and are finite, then we say pi is dominated by p2 in the Shannon transform sense, and 
write Pi <s P 2 to denote log(l -f p/s)pi(ds) < log(l -f p/s)g.2(ds), for all p > 0. 

Some properties of the Shannon transform order are the following: 

PI: If /“pi(ds) = /“p2(ds) < cx), then pi <Lt P 2 ^ Pi <s P 2 

P2: Pi <s P2 /o°° cr(-s)pi(ds) < o-(s)p2(ds), for all {a G TBX\a{ 0 ) = 0 }. 

P3: Let /“pi(ds) < /“p2(ds). Then pi <s p 2 a(s)pi(ds) < a(s)p2(ds), for 

all a G TBX. 

III. System Model 

The system model under eonsideration is the following: 

y = \/^b + w , (5) 

where y is the reeeived signal, p > 0 represents the average signal to noise ratio (SNR), b is the 
input symbol to the fading ehannel, w is eircularly symmetrie additive white Gaussian noise with 
zero mean and unit variance. In ([5]), X represents the magnitude square of the complex baseband 
equivalent fading coefficient, and is called as the instantaneous channel power. We normali z e 
E [X] = 1. Unifying the various distributions of X under the assumptions of multipath fading, 
shadow fading and composite multipath/shadowing is the focus of this work. The unification 
permits generalized performance analysis of the system with respect to metric such as the average 
SER, which is defined in Section IVI-AI 
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IV. Infinitely Divisible Fading Distributions 


A. Infinite Divisibility 

A probability distribution is said to be infinitely divisible (ID) if, for each n > 1, it can be 
decomposed into n identical convolution factors. Some examples of such distributions are the 
Gamma, the Poisson and the Gaussian distributions. In this paper, we refer to RVs with ID 
distributions as ID RVs. Non-negative ID RVs are characterized by the fact that the negative 
logarithm of the real-valued Laplace transform, i.e., the Laplace exponents of such RVs are 
Bernstein functions. The set of all nonnegative ID RVs is denoted by XV. In a wireless com¬ 
munications context, the resulting effect of coherent combination of contributions from multiple 
different sources can be viewed as fading nonnegative ID RVs. However, even shadow fading 
distributions like the log-normal which arise from infinite products of RVs turn out to be infinitely 
divisible. Indeed, almost all known instantaneous channel power distributions are ID, with 
examples listed in Section IIV-DI The infinite divisibility of the lognormal distribution, proved 
by Thorin [l^ revealed a very rich class of ID probability distributions known as generalized 
gamma convolutions (GGC), which possess extremely useful mathematical properties applicable 
to wireless systems. In addition, several popularly used fading distributions are not only ID, but 
also GGC. This fact motivates us to consider this class of ID distributions, which is discussed 
next. 

B. Generalized Gamma Convolutions 

The subset of non-negative ID distributions for which the Laplace exponent belongs to TBV, 
constitutes the set QQC of GGC probability distributions. In other words, the LT of the RV 
admits the form 



( 6 ) 


where the parameter a > 0, and the nonnegative measure /x(-) uniquely characterize the GGC 
distribution. In this paper, we call the Thorin measure of the Laplace exponent of the GGC (e.g. 
/i in ®) as the Thorin measure of the GGC. The degenerate GGC distribution is obtained with 
/i(dM) = 0. 





GGCs are closely linked with Gamma RVs having density 


fx{x) = 9^x^-^exp{-ex)/T{D), 


\D^D-l 


(7) 


for X E M"*", where 0 > 0 is known as the rate parameter, and D > 0 is known as the shape 
parameter, denoted as Gamma{D,9). In this context, a arbitrary GGC is a weak limit of a 
sequence of convolutions of gamma distributions with suitably chosen parameters. In wireless 
communications, the Thorin measure of a GGC has a physical interpretation: Consider the case 
where a continuum of fading channels with Gamma distributed instantaneous power is coherently 
summed up. If the sequence of shape parameters corresponding to each rate parameter m in [0, cx)) 
of the Gamma distributed channels is then the Thorin measure of the GGC channel is fi. 

In what follows, we will restrict our attention to nondegenerate GGCs with a = 0, because 
a merely represents a translation of the instantaneous channel power PDF, and performance 
analysis for such channels can be obtained from systems where a = 0 . 

The variation exponent of the LT of the GGC near infinity will correspond to the diversity 
order, and this quantity relates to the Thorin measure as follows: 

Lemma 1. /|^ Theorem 4.1.4] The LT of a non-degenerate GGC is regularly varying near oo, 
and the variation exponent equals the diversity order. That is, limp^oo f ifp) 14’ [p) = where 


OO 


D = //i(ds). 


0 


The variation exponent of the LT of a GGC near oo will have applications in asymptotic 
performance analysis, as discussed in Section |Vll The Thorin measure also determines the 
expected value of a GGC RV: 


Proposition 1. Let X be a GGC. Then E [X] = l/up{du), if the integral is finite. 


Proof: It is easy to show for any nonnegative RV that E[X] = — limp^o dfx (p) /dp. Since 
the LT of a GGC is of the form fix (p) = exp(— 7 (p)), E [X] = limp^.o d 7 (p)/dp, where 7 (-) 
is the Laplace exponent of the GGC. Representing 7 in terms of the Thorin measure. 


OO 


d/ log(l + s/u)p{du) 


E [Xl = lim -— 

p^o dp 


OO 



(8) 


0 




where /r is the Thorin measure. The limit and integral can be interchanged whenever u ^n{du) 
oo, because under this condition dominated convergence is satisfied. This completes the proof. 

■ 

If X represents the instantaneous channel power in a fading system, then Proposition [T] gives 
the average channel power of the system. As mentioned in Section Uni we assume that the mean 
of the GGC exists and is normalized to 1, as per the discussion after (|5]). In doing so, we ensure 
that p in ([5]) reflects the average signal to noise ratio of the system. 

It is straightforward to see that, if Xi and X 2 are two GGCs with Thorin measures pi and 
fi 2 , then Xi + A 2 is a GGC with a Thorin measure pi + p 2 - 

Some useful properties of non-degenerate GGCs are now presented. 

The density of a GGC always exists [| 8 l Theorem 4.1.2, Theorem 4.1.3], and can be written 
in terms of c.m. functions as follows: 

Lemma 2. /[^ Theorem 4.1.1, Theorem 4.1.2] The density of a GGC with diversity order D < 00 
can be written in the form 

fx (a:) = x^~^h{x) , (9) 

where h{-) is a c.m. function, which is slowly varying at 0. Moreover, h(0+) is finite if and only 
if /“logMp(dM) < 00 , and in this case 

logMp(du)^ . ( 10 ) 

The representation of the PDF of a GGC in the form (|9l) is particularly useful to obtain 
asymptotic performance metrics as shown in Section IVI-A[ GGCs are also mixtures of Gamma 
distributions, as summarized in the following Lemma: 

Lemma 3. Theorem 4.1.1] If X is a GGC with diversity order D < 00 , then X = AZ, 
where Z ~ Gamma{D, 1) and A > 0 are independent. Moreover, h{x) in ® is given by 

h{x) = {r{D))-^E[A-^exp{-x/A)] . (11) 

Interpreting X as the instantaneous channel power of a GGC channel, from Lemma [3] it is 
straightforward to see that every GGC channel with finite diversity order is equivalent to a 
variance mixture of a Nakagami-m channel. This is because the instantaneous channel power 


h(0+) = 


r(B) 


exp 
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corresponding to Nakagami-m fading is Gamma distributed If22l . It is shown in Seetion lTV-DI that 
several fading distributions sueh as Nakagami-m, Pareto and generalized gamma distributions 
have instantaneous ehannel powers whieh are GGC RVs, and the eorresponding diversity orders 
are finite. 

C. Hyperbolically Completely Monotone RVs 

A subset of GGC, eonsisting of the distributions for whieh the density is a hyperbolieally 
eompletely monotone (HCM) funetion is denoted by HCM (HCM c GGC C IV), and the 
RVs eorresponding to distributions in this set are ealled HCM RVs. A funetion g is defined to 
be a HCM funetion if g{uv)g{u/v) ean be written as (Sl p. 68] 

OO 

g{uv)g{u/v) = j ex.'p{—Xu{v + v~^))K{dX]u),K{dX]u)>0, (12) 

0 

for eaeh m > 0. HCM RVs enjoy the following properties under produets and quotients of RVs, 
whieh will be useful in obtaining elosure results for wireless systems involving multiple fading 
RVs: 

Lemma 4. If Xi, i = 1,2,... are independent HCM RVs, and Y is a GGC independent of Xi, 
then 

(i) Xi G HCM., for |g| > 1. 

(ii) X 1 X 2 G HCM. 

(in) XiY G GGC. 

(iv) e GGC, for \q,,f>l. 

D. Examples of ID Fading Distributions 

In this work, we view instantaneous ehannel power distributions as members of XV. To justify 
that the XV elass is indeed a unifieation of fading distributions, many of the popularly used 
fading ehannels are shown to have ID instantaneous ehannel powers. Among these ehannels, 
distributions for whieh the instantaneous ehannel power is a GGC or HCM RV are identified, 
as this faet will be useful in the analysis of systems, where the overall instantaneous ehannel 
power RV involves a sum or produet of two or more RVs (See Seetion IVlllh . Explieit examples 
of RVs whieh are ID, GGC and HCM are now listed. 






II 


1) Rician-K Fading: Rician-i^ distribution is used to model the fading envelope in certain 
line-of-sight scenarios ^l3[ p. 21]. For this case, if X is a nonnegative RV corresponding to the 
instantaneous channel power, then \/X is Rice distributed with parameter K. The LT of X is 
given by [|^ p. 19] 


(t>x (s) 


l + K ( Ks \ 


(13) 


It is now shown that X is an ID RV. To begin with, the Laplace exponent of X is obtained from 
(fT3]) as 

Using Frullani’s representation log(l + a;) = — exp(— xm)) exp{—u)/udu, for x > 0 lf24l 

p. 6], it is straightforward to show that 7x(s) can be written in the form ([T]) with a = b = 0 and 

r(s) = exp(-(l + K)s) + X) + ^^ ;. (15) 

Thus, 7 x(s) is a Bernstein function, and consequently, X is an ID RV. However, this is not a 
GGC because the Laplace exponent 7x(-) is not a Thorin-Bernstein function. To see this, recall 
from Section |n| that a Bernstein function with Levy measure r is a Thorin-Bernstein function, 
if and only if sr(s) is a c.m. function. However, by differentiating sr(s), where r(s) is defined 
in (fTS]) . it is seen that the first derivative is not non-positive at all values of s, which shows that 
st{s) is not c.m.. 

The instantaneous channel power corresponding to Rician-iF fading is therefore an ID RV, 
but not a GGC. 

2) Nakagami-m Fading: Nakagami-m fading envelop, which corresponds to a Gamma dis¬ 
tributed instantaneous SNR, is used to model line of sight channels. In addition, the Gamma 
distribution has also been used to model long-term fading effects in the literature [l25ll . Moreover, 
the exponential distribution, which is a special case of a Gamma RV with unit shape and rate 
parameters, arises as the instantaneous fading power distribution in Rayleigh fading channels. 

The instantaneous channel power X of Nakagami-m fading is Gamma distributed with shape 
m and rate m [|22l . In this case, X is a HCM, GGC and ID RV. To see that X is a GGC, observe 
that the LT of X is given by (px (s) = exp(—m log(H-s/m)), which satisfies the representation 
of the LT of a GGC specified in ®. It is then straightforward to show using ® that the density 
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of the Thorin measure in this ease, is given by 

dfi{u) = ml [{m}] , (16) 

and the Thorin mass is m, as ean be verified by rewriting the PDF in the form Further, h{-) 
is given by h{x) = mexp(—ma:)/r(m). The proof that X is a HCM RV is given in |[8l p. 75], 
and eonsequently a GGC and ID RV. 

3) Nakagami-q (Hoyt) Fading: Nakagami-g (Hoyt) fading is observed in seenarios sueh as 
satellite eommunieations subjeet to ionospherie seintillation Il26l . Speeial eases of Hoyt fading 
are Rayleigh fading (q = 1) and one-sided Gaussian (q = 0). 

The in-phase and quadrature eomponents of the eomplex baseband fading RV for the ease of 
Nakagami-g fading are independent zero mean Gaussian RVs with unequal varianees, sueh that 
the ratio of the quadrature varianee to the in-phase varianee equals g Il27l . Il2^ . Further, the 
average ehannel power is unity if the sum of the varianees is 1. As a result, the instantaneous 
ehannel power X ean be written as 


X = 


2 5 


(17) 


1 -f g2 1 -f g2 

where Wi ~ A/'(0,1) and kF 2 AA( 0 ,1) are independent, and 0 < g < 1. 

It is now argued that X is a GGC and henee also ID. Toward this end, first observe that 
Xj, i = 1,2 is a GGC, beeause its Laplaee exponent is given by 7 ( 5 ) = (l/2)log(l -f 2s), 
whieh is a Thorin-Bemstein funetion. It then follows from (fTTI) that X is a GGC, beeause GGCs 
are elosed under nonnegative sealing of RVs and addition of RVs. Moreover, using the identity 
4>ax (s) = 4>x (as) for any nonnegative RV X, it ean be shown that the density of the Thorin 
measure is given by 


dqb(u) 1 

du 2(1+ q‘^)' 


d-^} 


+ 


T 


2(1+ g2)' 


{i±^} 

^ 2g2 \ 


(18) 


and the Thorin mass is 1. It is also noted that X is not HCM, beeause the Thorin measure of 
a HCM RV is allowed to have at most one atom [[8l p. 88], whieh is not the ease for X, as 
observed from (fT^ . Therefore, X for Hoyt fading is a GGC and ID RV, but not a HCM RV. 

4) General Lognormal Distribution: The lognormal distribution is eommonly used to model 
the long-term shadowing effeet in wireless ehannels [|2^ . 

The density of the general lognormal RV is given by 

f (loga:)^' 




exp 


2cr2 


fx (a;) = cx‘ 


(19) 
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with yU G M, cr > 0, and c > 0 is a normalizing constant. It has been shown in [|8l p. 74] that the 
general lognormal distribution belongs to 1-LCM., and hence also QQC and XV. Although the 
Thorin measure is unknown, and the diversity order can be calculated to be cxd (H p. 74]. 

In wireless systems, when multiple channels with the lognormal shadowing and XLCM. fading 
effects combining together on each individual channel, the total received channel gain using 
MRC can be shown to be GGC RV using Lemma |4| 

5) Generalized Gamma (Stacy) Distribution: The instantaneous channel power X is general¬ 
ized Gamma distributed, when the envelope of the fading amplitude is modelled as a generalized 
Nakagami-m RV [l25l . The Weibull distribution is a special case of the generalized Gamma 
distribution, and has been used to approximate the multipath wireless channel from channel 
measurements [|29l . Other distributions which are special cases of the generalized Gamma 
distribution are the inverse-Gaussian distribution and the Nakagami-m distribution. Moreover, 
the generalized gamma distribution arises as the instantaneous channel power distribution in 
Nakagami-m fading channels with receiver non-linearity, where the non-linearity is captured as 
a power parameter, since X can be written as ^ p. 13] 

a: = , (20) 

where Y is Gamma distributed with shape parameter e/r and rate C2, denoted as V ~ Gamma(e/r, C2). 

The density of X is given by 

fx (x) = ciy ^ exp (-C2X'’) , (21) 

where 0 < |r| < 1, ci, C2 > 0, and e/r > 0. 

It is now shown that this distribution belongs to XiCM. and consequently QQC and XV. Recall 
that the Gamma distribution belongs to XiCAi, as shown in Section ITV-D2I Using Lemma HI it 
is then seen that the generalized gamma distribution also belongs to XLCAi, and is therefore a 
GGC. While the Thorin measure for the general case is unknown, the diversity order is equal to 
e if r > 0, and 00 otherwise, as can be verified by representing (|2TI) in the form (|9|), and further, 
the function h{-) in (Hj) is given by h{x) = ciexp{—C 2 X^). 

6) Product of Generalized Gamma Random Variables: The product of independent general¬ 
ized Gamma RVs has been proposed as a unified fading model in [l30l . As a special case of this 
family of RVs, product of N independent Gamma RVs is obtained, which is the instantaneous 
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channel power distribution in the generalized-it' fading model OTl . Another speeial ease of this 
distribution is the A^^Nakagami fading model [O, [|32l, where the envelope is a produet of N 
independent but not neeessarily identieally distributed Nakagami-m RVs. 

Consider X to be a product of N independent but not necessarily identieally distributed 
generalized Gamma distributed RVs. This distribution is HCM as observed from an applieation 
of Lemma |4] and the fact that the generalized Gamma distribution belongs to T-LCAi. As a 
eonsequenee, the distribution is also in QQC and XT). 

7) Spherically Invariant Random Process: The spherieally invariant random proeess (SIRP) 
model has been proposed as a unified model for the fading envelope distribution in the literature 
ifTl . In this ease, the instantaneous ehannel power X is given by 


X = AE , 


( 22 ) 


where E is an exponential RV, and A is a positive valued RV. 

It is now shown that X for this case is always an ID RV, and is GGC and HCM for certain 
speeial cases. To see that X is ID, observe that X is a varianee mixture of exponentials, 
whieh is ID, aeeording to dH Theorem 2.4.3]. Not every member of SIRP is a GGC. A simple 
eounterexample is the ease of Rician fading, which is a member of SIRP according to [[II, and 
not a GGC as seen in Seetion ITV-Dll The speeial eases of SIRP whieh are GGC RVs and HCM 
RVs are now eonsidered. X is a GGC, whenever A is a GGC, aeeording to Lemma |4| Also, X 
is a HCM RV for all such members of SIRP where A is a HCM RV. 

8) Positive Stable Distribution: Positive stable distributions have been used to model the 
interferenee at the primary reeeiver in a cognitive radio network, when the interfering seeondary 
terminals are distributed in a Poisson field, and there is a guard zone around the primary reeeiver 

m. 

The positively skewed stable distribution is a heavy tailed distribution, whieh is eharaeterized 
by its LT as (j)x (s) = exp(—s^), 0 < r < 1. No elosed form expression for the distribution or 
density is known in general. Nevertheless, it is known that this distribution is a GGC dH p. 35], 
and the density of the Thorin measure is given by 


dp,{u) 


r sill rvr 


111 / /I 1 

- u 

TV 


(23) 


dti 


and diversity order equal to oo. 
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9) Pareto Distribution: Pareto RVs are heavy-tailed distributions used to model signal to 
interferenee ratios in interferenee dominated seenarios. A Pareto RV ean be written in the form 


m p. 14] 



(24) 


where Xj ~ Gamma{kj, 1), j = 1,2 are independent. The faet that Pareto RVs are ratios of 
Gamma RVs raised to a non-negative power, as shown in (l24l) leads to many wireless systems with 
Pareto instantaneous ehannel power. For example, eonsider a system with a transmitter, reeeiver 
and an interfering terminal. Suppose the ehannel between the transmitter and the reeeiver, and 
that between the interferer and reeeiver are both Nakagami-m ehannels (with possibly different 
parameters), and the reeeiver is non-linear with a power nonlinearity. Then from (l24l) . it is elear 
that the instantaneous ehannel power is a Pareto RV. It is also noted that speeial eases of the 
Pareto distribution with (ki = k 2 = l,r > 0) have been used to model the instantaneous signal 
to interferenee power in interferenee dominated networks Q. 

The density of a Pareto RV is given by 


(25) 


B{k2,ki) (1 + a:^)^i+'=2 ’ 


where B{a, h) := is the Beta funetion, r e [—1,1], kj > 0, j = 1, 2. 

The Pareto distribution belongs to PLCAi. This is beeause Gamma RVs are HCM, and using 
Lemma |4| on (|2^ . it follows that a Pareto RV is HCM. The Thorin measure is unknown in the 
general ease, however the diversity order is obtained as kir [[8l p. 74]. 

As diseussed above, almost all of the fading distributions with ID instantaneous ehannel power 
also belong to QQC. It is therefore reasonable to foeus on ehannels with instantaneous ehannel 
power in QQC. Hereafter, we refer to ehannels with GGC distributed ehannel power as GGC 
ehannels. 


V. GGC Fading Channels as Gamma Mixtures 


Let X be a GGC with finite diversity order D. As diseussed in Lemma [U X ean be written 
as AZ, where Z is Gamma distributed with parameters (77,1), and A is a nonnegative RV 
independent of Z. Using the property tZ ~ Gamma{D,l/'t) applieable to Gamma RVs, it 
is easy to see that X ean be written as AZ, where Z is Gamma distributed with parameters 
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{D,D), and yl is a nonnegative RV independent of Z. This representation shows that a GGC 
with ehannel with finite diversity order is a Nakagami-m fading ehannel with shadowing, where 
m = D. This is beeause the instantaneous ehannel power of a Nakagami-m ehannel with m = D 
has the same distribution as Z. Moreover, the average ehannel power of the shadowing RV is 
1, sinee it has been assumed that E [X] = 1. 

The distribution of the shadowing RV A is not trivial to obtain, when only the distribution of X 
is known. In this ease, we adapt the approaeh used to obtain the mixing distribution for mixture 
of exponential distributions [|T1|. This method is based on the faet that the Mellin transform of 
a produet of independent RVs is the produet of the Mellin transforms of the RVs, and that the 
distribution and its Mellin transform pair are unique. The Mellin transform of a nonnegative RV 
A is defined as 


^ a ( s ) = j , (26) 

0 

and M^^{s) = M^(s)M^(s). Therefore, the PDF of A ean be obtained using the inverse Mellin 
transform as follows: 


f^(x) = M 


-1 


Mx(s) 


(27) 


The Mellin transform of Z ean be evaluated as m^“^r(s + D — 1)/T{D) [l34l p. 312]. Thus, 
substituting the Mellin transform of X in (1271) . the density of A ean be obtained. Referenee [[U 
uses this approaeh to obtain in terms of Meijer-G funetions . 


VI. Performance Metrics for GGC Channels 

In [l35l . the diversity and the asymptotie error rate of a wide range of ehannel distributions 
are analyzed using the theory of regular variation. It has been proved that the diversity order is 
equivalent to the variation exponent of the CDF of the ehannel power gain at the origin. The 
asymptotie average error rate eharaeterization applies to a general set of ehannel distributions. 
In this Seetion, with the establishment of the infinitely divisible RVs, asymptotie expressions for 
average SER of GGC ehannels and/or the funetion h{-) defined in dH), and the diversity orders 
in terms of the Thorin measure are obtained. 
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A. Asymptotic Symbol Error Rate 

The average SER (p) is defined as follows: 

P. (p) = E |P,(pX)] , (28) 

where Pe{pX) is the instantaneous SER, and is dependent on the eonstellation of ehoiee. 

The analysis of the asymptotie SER of M-PSK and DPSK are eonsidered. In the asymptotieally 
high SNR regime, Pg (p) is given by 

Pe (P) ~ CQh(p-^)p-^ ,p^oo , (29) 

where h{x) is slowly varying at 0 ll^ . Using (l29l) . generalized definitions quantities related to 
asymptotic SER are obtained. The diversity order D is the asymptotic slope of the average SER 
on a log-log plot versus the average SNR, and is defined as the variation exponent of (1281) at 
oo [1351. The term cq in (l29l) is a modulation dependent constant, and is given by cg = vr/d for 
the case of DPSK |[35l . and 

^{1—iW~^)7r 

Cq = (7rsin^'^(M“V))“^ / sin^^ 6*d6', (30) 

Jo 

for the case of MPSK lf35l . 

Now the average SER performance of a GGC channel at high SNR is now considered, in 
order to obtain the diversity order in terms of the Thorin measure. 

Proposition 2. Let a GGC fading channel have diversity order D < oo. Then for DPSK/MPSK, 
we have Pg (p) ~ CQh{p~^)p~^, as p ^ oo, where cq depends on the modulation scheme, and 
hf) is as defined in ® or (ITTI) . 

Proof: According to [[35l Theorem 4], Pg (p) can be written as CQh{p~^)p~^, as p ^ cx) for 
DPSK and MPSK, if the asymptotic density of X has the form fx (x) = x^~^h{x), as x —)• 0. 
GGCs with finite diversity order satisfy this condition as seen from Eemma[2l ■ 

It is observed from Proposition [2| that the diversity order of a GGC channel is as follows. 

Corollary 1. The diversity order of DPSK/MPSK over a GGC fading channel with diversity 
order D < oo is equal to D. 

The asymptotic SER obtained in Proposition [2| applies to any GGC fading channel with finite 
diversity order, and shows the dependence of the parameters of the GGC on the asymptotic 
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SER parameters. For GGC ehannels satisfying a eertain condition on the Thorin measure, the 
asymptotic SER can be written in the form (p) ~ Gp~^, where G is known as the array 
gain, which represents the shift of the average SER curve to the right on a log-log plot 


Proposition 3. Let the Thorin measure of a GGC fading channel with finite diversity order D 
satisfy /j°°(logM)/i(dM) < oo. Then (p) ~ {Gp)~^, as p ^ oo, where 


G = 


Proof: From Proposition |2l 


r(Z)) 


exp 


{\ogu)p{du) 


■ D 


(31) 


CopPe(P) 
lim ^ ^ = 1 . 


(32) 


The limit of the ratio in 


p-i-oo h(p 

can be written as a ratio of limits of the numerator and denominator. 


if and only if J“(logM)/i(dM) < oo, because this is the condition under which limp_j.oo h{p 
oo, as can be seen through an application of Eemma [2l Therefore, (|3^ simplifies to 


-l^ 


< 


lim (P) = lim h{p ^) . 


(33) 


The right hand side of (l3^ is given by T{D) ^ exp{J^{logu)p{du)), as seen from EemmalU 
Substituting this in (l3^ and rearranging, 

PeiP) 


lim 


= 1 . 


(34) 


p-^oo p ^exp(f^ (logu)p(du)) 

This proves the Proposition. ■ 

According to Proposition [3l the slow varying function h(x) in Proposition[2]becomes a constant 
at 0. This conclusion is not necessarily true for fading distributions which are not GGCs, or are 
GGCs which do not satisfy /“(logM)/i(dM) < oo. An example of a fading distribution for which 
the slow varying function h does not become a constant at 0 is the case of generalized iT-fading, 
as discussed in [|35]|. 

It can be seen from (l3TI) that the array gain is inversely proportional to /i(0+), where h(-) 
is as defined in dH), because the term r(ii))“^ exp(/Q°°(logM)p(dM)) equals h(0+), according to 
Eemma [2l The array gain will be useful in calculating the performance difference between two 
GGC fading channels at high SNR, as discussed in Section IVI-BI 
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B. Quantifying Asymptotic SER Performance Gaps 

In the absence of a unified fading model, it is difficult to obtain the SNR gain of one fading 
distribution with respect to another, without explicitly obtaining the asymptotic SER expressions 
for both channels. In what follows, the properties of GGCs make it possible to produce a closed 
form expression for the SNR gain of one GGC with respect to another, when comparing the 
SERs at high SNR. While stochastic ordering techniques, which will be illustrated in Section 
IVIIl are capable of yielding GGC channel comparisons based on very general performance 
metrics without the need for closed form expressions at every value of average SNR, an inherent 
limitation of this method is that it is not possible to quantify the exact performance gap of one 
channel with respect to the other. 

The SNR gain of one channel with respect to another is defined as the difference between 
the SNR (in dB) for the asymptotic SER of the two channels to be equal, and is denoted as 
QdB (-^ 1 ,-^ 2 )- Mathematically, 



(35) 


Pe ^ (P) 


where (■) and P^^^ (■) are the SERs of the channels with instantaneous channel power Xi 


and X 2 respectively. 

Proposition 4. Let Xi and X 2 correspond to the instantaneous channel powers of two GGC 
channels satisfying \ogupj{du) < 00 , and dj(du) = D, j = 1, 2. If MPSK or MQAM is 
employed, the SNR gain is given by 


00 



(36) 


0 

where Pj{-) is the Thorin measure of Xj, j = 1,2. Equivalently, 



(37) 


where hj{-) are obtained from the canonical GGC PDE ® of Xj,j = 1,2. 


Proof: Relation (1371) is proved first. Then (l3^ is obtained using the relation (fT^ . Eet pi 
and p 2 denote the average SNRs for channels Xi and X 2 required to obtain the same asymptotic 
SER. In other words, {Gipi)~^+0{pi^~^) = {G 2 P 2 )~^+0{pf^~^). Relation (1371) then follows 
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by ignoring the higher order terms, taking logarithms both sides, and substituting the array gain 
expressions from (ISTl) . ■ 

VII. Stochastic Ordering of GGC Distributions 
A. Laplace Transform Ordering of GGC 

The performanee of wireless systems are quantified by averaging a metrie (e.g. bit or symbol 
error rates, or ehannel eapaeity) over the distribution of the random ehannel. Quantifying the sys¬ 
tem performance relies on single parameter comparisons between channels using characteristics 
such as diversity order [|^ . This approach does not provide a unified frame work to compare 
channels across many different performance metrics. The theory of stochastic orders provides 
a comprehensive framework to compare two RVs. It allows comparing systems using different 
metrics with properties such as monotonicity, convexity, and complete monotonicity, which shed 
light into the connections between performance metrics such as error rates and ergodic capacity. 

In this Section, the stochastic ordering of GGC fading channels is considered. This will help 
in comparing two GGC fading channels based on general performance metrics which are either 
completely monotone (such as SERs), or possess a completely monotone derivative (such as the 
ergodic capacity). Toward this end, the Laplace transform ordering framework proposed in ll3^ 
is employed. LT ordering between a pair of instantaneous channel power distributions implies 
that the average SER of a constellation with c.m. SER, such as the case with MPSK and MQAM, 
will be ordered at all values of average SNR. 

In the context of LT ordering of GGCs, a duality between the Shannon transform ordering of 
the Thorin measures with the LT ordering of the fading distributions is straightforward to see: 

Proposition 5. Let X and Y be two GGCs with Thorin measures px and py- We have X <Lt 
V px <s Ty- 

Proof: The proof of this proposition follows from (01) applied to GGCs. ■ 

The connection between the LT order and the Shannon transform order as suggested by Propo¬ 
sition [5] can be exploited in obtaining new ordering relations between GGC fading distributions, 
by using the properties of one stochastic order to benefit the other. Lor instance, by observing that 
g{x) = log(l -f px) is a Bernstein function for p > 0, it is seen that px <Lt Py dx <s Ty, 
for cases when the Thorin masses of X and Y are identical. Now, using Proposition [51 it is 
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concluded that nx <Lt X <Lt Y, if the diversity orders of X and Y are equal. 

Therefore, the generalized LT ordering of the Thorin measures of two GGCs with equal diversity 
order implies that the average SER performance at all SNR for the first one is better than that 
of the second at all SNR. 

Laplace transform ordering of GGCs can also be obtained through the observation that any 
GGC with finite Thorin mass is a gamma variance mixture. This leads to the following Propo¬ 
sition: 

Proposition 6. Let X = AUi and Y = BU 2 be two GGCs, where A and B are nonnegative 
RVs, independent of Uj ~ Gamma{mj,l), j = 1,2. Then mi < m 2 and A <Lt B implies 

X <Lt 

Proof: To begin with, observe that the LT of a Gamma RV monotonically decreases with 
the scale parameter, when the rate is unity. Therefore, if Ui and U 2 are Gamma RVs with shape 
mi and m 2 respectively, and unit rate, then Ui <Lt U 2 rni < m 2 . Let A and B be 

non-negative RVs independent of f/i and 1 / 2 , which satisfy A <Lt B. Then 

E [exp{-pAUi)] = Eyi [Et/j [exp{-pAUi)]] (38) 

> Ea [Eu 2 [exp{-pAU 2 )]] (39) 

> Eb [E [/2 [exp{-pAU 2 )]] , (40) 

where (l39l) is obtained from (1381) by observing that Ui <Lt U 2 , and exp(—x) in (l3^ is a c.m. 
function. Similarly, (l40l) follows from (l39l) . as A <Lt B. ■ 

While Proposition [5] connects the ordering of the Thorin measures and the LT ordering of the 
corresponding GGCs, Proposition enables the LT ordering of pairs of GGC RVs for which the 
Thorin measure may not be available in closed form. 

B. Comparison of GGC Channels with Equal Diversity Orders 

It is well known that the AWGN (no fading) channel is a benchmark to the performance of any 
fading channel, with respect to symbol error rates and ergodic capacity. Intuitively, this is because 
the AWGN channel has infinite diversity order (since a Nakagami-m with m —)■ 00 is an AWGN 
channel). However, when two GGC fading distributions with the same diversity order are to be 
compared, it is in fact observed that the Nakagami-m fading scenario is the best possible fading 


22 


channel with respect to symbol error rates of 1-dimensional or 2-dimensional constellations 
(whieh are eonvex funetions) and the ergodie eapaeity (whieh is a eoneave funetion), as seen 
from the following Proposition. 

Proposition 7. Let X be a GGC with Thorin mass D < oo and EX = 1. Then U <cx X, where 
U ^ Gamma{D^D). 

Proof: Aeeording to Lemma[3l X ean be written as X = AU, where U is Gamma distributed 
with shape D and rate D, and A is a nonnegative RV independent of U. Now, let g he a eonvex 
funetion. Then Ef/ [5f(E [A] U)] < E [g{X)], as a eonsequenee of Jensen’s inequality. E [A] = 1, 
beeause E [X] = 1, by assumption. The proof is thus eoneluded. ■ 

It is now possible to see from Proposition |7] that among all GGC fading ehannels with 
a given diversity order and unit average power, the Nakagami-m fading ehannel forms the 
benehmark ehannel with respeet to eonvex performanee metries sueh as SERs of 1-dimensional 
or 2-dimensional eonstellations, and eoneave metries sueh as the ergodie eapaeity. This is beeause 
the instantaneous ehannel power of a Nakagami-m ehannel is Gamma distributed with shape m 
and rate m. In other words, if the diversity order of a GGC is D, then the average SER of any 
1-dimensional or 2-dimensional eonstellation is lower bounded by that of a Nakagami-m ehannel 
with m = D, and further, the ergodie eapaeity is upper bounded by that of the Nakagami-m 
ehannel with m = D. 


VIIE Systems Involving Multiple GGC RVs 

Many systems of praetieal interest involve eombinations of multiple fading eomponents. In 
what follows, the conditions under whieh the overall end-to-end instantaneous ehannel power is 
a GGC with finite diversity order is obtained. This permits the eomparison of two GGC ehannel 
systems using Proposition 0 This faeilitates the use of Proposition |4] to quantify the high SNR 
gain of one GGC ehannel with respeet to another, sinee this propositions require finite diversity 
order. 

A. Composite Fading Systems 

Many fading models in wireless eommunications attempt to capture both the short term fading 
effeet and the long term shadowing effeet through a eomposite fading distribution. In sueh oases. 
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the instantaneous ehannel power is modelled as a produet of two RVs, one corresponding to the 
short-term effect, and the other corresponding to the long-term fading effect. In other words, 
the overall instantaneous channel power Z can be written as Z = XY. For such systems, Z 
is a GGC if X is a GGC and Y is HCM, as seen through an application of Lemma |4l It can 
be observed that the diversity order of Z is finite if the diversity order of either X or L is 
finite. As a result, common composite models such as the Rayleigh-lognormal model [l37]| . the 
Nakagami-lognormal model [l38l . Weibull-Gamma composite model [|39l, Generalized-X [|3TI . 
and the Gamma-shadowed generalized Nakagami fading model [l25l have GGC instantaneous 
channel powers with finite diversity order. 

It can then be concluded using Proposition that, for two different composite fading distri¬ 
butions Z\ and Z 2 , each with finite diversity order, if the diversity order of Z\ is larger than 
that of Z 2 , then to establish LT ordering of the composite distribution, it is sufficient to establish 
LT ordering of the mixing distributions. Further, if the Thorin measures of the two composite 
distributions satisfy log (dw) < 00 , then Proposition |4] can be used to quantify the high 
SNR gain of one system versus the other. 

B. Diversity Combining Systems 

Consider a single-input multiple-output diversity combining system with L receive antennas, 
and complete CSI at the receiver. It is now proved that under different assumptions on the 
instantaneous channel power of each branch, the end-to-end instantaneous channel power of 
combining scheme such as MRC. In the case of MRC, the instantaneous end-to-end channel 
power Xmrc = J2i=i ^ GGC whenever Xi is a GGC, since a sum of GGC RVs is a 
GGC. It is therefore straightforward to see that the average SER of DPSK is given by one-half 
0 (p) in ®, where the Thorin measure is the sum of Thorin measures of each component of 
the MRC system. The ergodic capacity and average SERs of other 2-dimensional modulations 
may however not be tractable in general, for example when the branches are independent but 
not identically distributed GGCs. In this case, it is possible to obtain a performance lower bound 
if the diversity order is known. To this end, if every component has finite diversity order, then 
the diversity order of the MRC system is finite, and is given hy D = I Now, 

since Xmrc is a GGC with diversity order D < 00 , using Proposition |71 it is seen that the 
performance of the MRC system is lower bounded by that of a Nakagami-m channel with 
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m = D, with respect to all convex metrics such as average SERs of 2-dimensional modulations 
and concave metrics such as the ergodic capacity. 

IX. Simulations 

In this section some of the theoretical results are corroborated using Monte-Carlo simulations. 
In Figure [H the performance of a 2-branch MRC system where the instantaneous channel powers 
are Pareto distributed with parameters (1,1,1.25) is simulated, and compared with that of a 
SISO Nakagami-m channel with m = 2.5. The average powers of both the systems have been 
normalized to unity. This simulation demonstrates that the average SER of the MRC system 
with Pareto distributed branches, which are GGC’s with finite diversity order is lower bounded 
by that of a Nakagami-m channel, as suggested in Section IVTTT-Rl 

In Figure [2l the performance of a Nakagami-m channel with m = 2 is compared with that of a 
fading channel with generalized gamma distributed instantaneous channel power with parameters 
(e = 2, r = 2, C 2 = 1). The parameters of the two distributions have been chosen such that the 
diversity order is 2 for both cases. The high SNR gain in dB obtained from the simulation is 
found to be ~ 1.7 dB. This agrees with the theoretically suggested value of 1.505 dB obtained 
from (1^ . with /ii(O-l-) = 4 and h2{0+) = 2. 

A simulation to provide an intuitive understanding of the structure of the Thorin measure 
and its effect on the average SER performance has been provided in Figure |4l To this end, two 
different GGCs X = Xi + X 2 and Y = Yi + Y 2 , where Xi,X 2 are independent and gamma 
distributed with parameters (2,1) and (2, 2) respectively, and Yi, Y 2 are independent and gamma 
distributed with parameters (1,1/2) and (3, 3) respectively, are chosen. In this case, the densities 
of the Thorin measure are as depicted in Figure [H The average SER for X and Y are obtained 
in figure IH It is observed that the average SER of X is consistently less than that of Y at all 
p > 0. Therefore, it can be inferred that the smaller the support of the density of the Thorin 
measure, the better the average SER performance is. 

X. Conclusions 

In this paper, it has been shown that the class of non-negative ID RVs serves as a unification 
of commonly used fading distributions such as Rician-7f, Nakagami-m and Weibull distribution. 
Furthermore, several unified fading distributions proposed in the literature, such as the Stacy 
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Fig. 1. Average SER performance of 2-branch MRC with Pareto(1.25) distributed SNR compared with that of Nakagami-m 
SISO channel with m — 2.5 



Fig. 2. Comparison of DPSK performance of Nakagami-m (m = 2) channel and generalized gamma distributed channel 
Gamma(l, 1)^^^ with equal avg. channel powers and diversity orders. 


distribution and the SIRP are shown to be ineluded in the elass of ID RVs. In addition, it is 
shown that several ID fading distributions are also GGCs. The properties of ID and GGC RVs 
find applieations in the performanee analysis and eomparisons of wireless systems. Speeifieally, 
the asymptotie average SER is obtained in a eanonieal form, and its depenenee on the Thorin 
measure of the GGC is revealed. On the other hand, eomparing two GGC ehannels in the LT 
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Fig. 3. Density of Thorin measures for two GGCs. X =Gamma(2,1)+Gamma(2, 2); Y =Gamma(l, 0.5)+Gamma(3, 3). 
Diversity order = 4, E [X] = E [F] = 3 



Fig. 4. SER performance of DPSK over two GGCs which are sums of 2 gamma RVs with different supports of Thorin 
measures. The less spread apart is, the better the performance. 


order is shown to be simplified, by revealing the eonneetion between the Shannon transform 
ordering of the Thorin measures and the LT ordering of the GGCs. This faeilitates eomparisons 
of the average SER or ergodie eapaeity of the two GGC ehannels. The properties of GGC RVs 
also permit us to quantify the asymptotie SER gap for any modulation seheme. Applieations of 
properties of GGCs in diversity eombining systems are also provided. 
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